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Abstract 

Protons and antiprotons at collider energies are a source of high energy Weizsacker- Williams 
photons. This may open a possibility to study exclusive photoproduction of heavy vector mesons at 
energies much larger than possible at the HERA accelerator. Here we present a detailed investiga- 
tion of the exclusive J/ip photoproduction in proton-proton (RHIC, LHC) and proton-antiproton 
(Tevatron) collisions. We calculate several differential distributions in ti,t2,y,4>, as well as trans- 
verse momentum distributions of J/^f's. We discuss correlations in the azimuthal angle between 
outgoing protons or proton and antiproton as well as in the (ti , ^2) space. Differently from elec- 
troproduction experiments, here both colliding beam particles can be a source of photons, and we 
find large interference terms in azimuthal angle distributions in a broad range of rapidities of the 
produced meson. We also include the spin-flip parts in the electromagnetic vertices. We discuss 
the effect of absorptive corrections on various distributions. Interestingly, absorption corrections 
induce a charge asymmetry in rapidity distributions, and are larger for pp reactions than for the pp 
case. The reaction considered here constitutes an important nonreduceable background in recently 
proposed searches for odderon exchange. 
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I. INTRODUCTION 



The diffractive photoproduction of J/^-mesons has been recently a subject of thorough 
studies at HERA [J @], and serves to elucidate the physics of the QCD pomeron and/or 



the small-x gluon density in the proton (for a recent review and references, see |3|). Being 
charged particles, protons/antiprotons available at e.g. RHIC, Tevatron and LHC are a 
source of high energy Weizsacker-Williams photons, and photoproduction processes are also 
accessible in hadronic collisions. Hadronic exclusive production mechanisms of mesons at 
central rapidities in pp collisions were intensively studied in the 1990-ies at energies of a few 
tens of Ge V 0] , and raised much theoretical interest for their potential of investigating exotic 
hadronic states (see e.g. [5|). Recently there was interest in describing diffractive exclusive 
production of heavy scalar jg, 0] and pseudoscalar [s[ mesons in terms of off-diagonal uninte- 
grated gluon distributions, which m ay p rovide insight into the related diffractive production 
mechanism of the Higgs boson (0, llCl| | and references therein). A purely hadronic mecha- 
nism for the exclusive production of J/^f mesons in proton-proton and proton-antiproton 
collisions was suggested as a candidate in searches for yet another exotic object of QCD, 
the elusive odderon exchange [III 13] ■ In order to identify the odderon exchange one has to 
consider all other possible processes leading to the same final channel which in the context of 
the searches for the odderon will constitute the unwanted background. One of such processes 
(and perhaps the only one at the level of fully exclusive J/\l/-production) is pomeron-photon 



or photon-pomeron fusion [12], [13|, [14|, which we study in this communication at a more 
detailed level than available in the literature. We feel that its role as a background for odd- 
eron searches warrants a more detailed analysis including energy dependence and differential 
distributions of the photoproduction mechanism in hadronic collisions. As will be discussed, 
the process considered here is interesting also in its own rights. 

An important concern of our work are absorption effects. More often than not absorption 
effects are either completely ignored or included as a multiplicative reduction factor, which 
is simply wrong for many observables (like distributions in t\ or £2) as we shall show in our 
paper. We think this point requires broader public spread as it often appears to be forgotten 
or ignored. We present a detailed analysis of several differential distributions in order to 
identify the absorption effects. We also put a special emphasis on interference phenomena. 
We will discuss also more subtle phenomena like the spin flip in the electromagnetic vertices 
and a charge asymmetry by comparing differential distributions in proton-proton and proton- 
antiproton exclusive J/ip production. 

In this work, we do not include a possible odderon contribution. We wish to stress, that 
the photoproduction mechanism of exclusive J/^'s must exist without doubt, and does not 
die out as energy increases. A related purely hadronic (Odderon) contribution with the 
same properties has not been unambiguously identified in other experiments, and hence 
cannot be estimated in a model-independent way. While certain QCD-inspired toy-models 
for CP-odd multigluon t-channel exchanges exist, they do not allow reliable calculations 
of hadronic amplitudes. In practice the magnitude of the corresponding Born amplitude 
strongly depend on the details of how to treat gluons in the nonperturbative domain, as 
well as on the modeling of proton structure. Furthermore, the energy dependence of the full 
(beyond the Born approximation, and beyond perturbation theory) amplitude is unknown 
and it can even not be excluded that this contribution would vanish with rising energy. It is 
not the issue of our paper to further discuss such models, we rather think that in the search 
for an odderon one should take further initiative, if substantial deviations from the more 
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conservative physics discussed here are found. 



II. AMPLITUDES AND CROSS SECTIONS 
A. 2^3 amplitude 



Here we present the necessary formalism for the calculation of amplitudes and cross- 
sections. The basic mechanisms are shown in Fig{TJ 




FIG. 1: The sketch of the two mechanisms considered in the present paper: 7IP (left) and IP 7 
(right). Some kinematical variables are shown in addition. 

The distinctive feature, when compared to photoproduction in lepton-hadron collisions, 
is that now both participating hadrons can serve as the source of the photon, and it is nec- 
essary to take account of the interference between the two amplitudes. Due to the Coulomb 
singularities in the photon exchange parts of the amplitude, the electromagnetic vertices 
involve only very small, predominantly transverse momentum transfers. Their effect is fully 
quantified by the well-known electromagnetic Dirac- and Pauli form factors of the nucleon. 
Regarding the photoproduction amplitude, we try to be as far as possible model indepen- 
dent, and take advantage of the precise knowledge of diffract ive vector meson production 
over a broad energy range available from experiments at HERA. The amplitude for the 
2 — > 3 process of Figfj] can be decomposed as 

Mllh 2 2hlhlv V ( s ' s i> s 2, h, t 2 ) = M lJP + M Wl 

= Wv K\J»\Pn \iK( qi , Xv)^^M% X 2 2 S v X h v 2 X \s 2 , t 2 , Q\) 

K\Mp2, A 2 ) e ;(g 2 , Xv)^^M x ;:l x i !T v x :;\ Sl , h, Ql) . 

(2.1) 

The outgoing protons lose only tiny fractions Z\,z 2 <C 1 of their longitudinal momenta. In 
terms of their transverse momenta p 1 2 the relevant four-momentum transfers squared are 
ti = —{pi + zfm^)/{\ — Zi) ,i = 1,2, and S\ ~ (1 — z 2 )s and s 2 ~ (1 — Zx)s are the 
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familiar Mandelstam variables for the appropriate subsystems. Due to the smallness of the 
photon virtualities, denoted by Q 2 = —ti, 1 it is justified to neglect the contribution from 
longitudinal photons, recall that (Tl{i*P ~^ J /^p) / &t{i*P — > J/^p) oc Q 2 /m 2 ,^ j^, Ill- 
Then, the amplitude for emission of a photon of transverse polarization Ay, and transverse 
momentum q x = —Pi, entering eq. fl2.ll) reads: 



(2.2) 



here = —(Xe x + ie y )/y/2, n\\e z denotes the collision axis, and <ti/2 is the spin operator 
for nucleon 1, \\ is its spinor. Fi and F 2 are the Dirac- and Pauli electromagnetic form 
factors, respectively. Here we have given only that part of the current, which gives rise to the 
logarithmic dz/z longitudinal momentum spectrum of photons, which dominates in the high- 
energy kinematics considered here. It is worthwhile to recall, that for a massive fermion, 
that includes a spin-flip contribution originating from its anomalous magnetic moment, 
k p = 1.79. Notice its suppression at small transverse momenta. The parametrization of 
the photoproduction amplitude which we used in practical calculations can be found in the 
Appendix. Above we already used the assumption of s-channel-helicity conservation in the 
7* — > J/ty transition, which for heavy vector mesons is indeed well justified by experiment 2 
In summary we present the 2 — ^ 3 amplitude in the form of a 2-dimensional vector 

as 

2 p 

M (Pi,P2) = ei — -rF\'\ 1 {Pn t i)M rh2 ^ Vh2 {s 2 ,t 2 ,Ql) 
2 p 

+e2—-^F\i\ 2 {P2,t 2 )M T h 1 ^vh 1 { s ii t iiQl) • ( 2 - 3 ) 

The differential cross section of interest is given in terms of M as 

da = 51 ^ 4g2 |M| 2 dydt x dt 2 d4> , (2.4) 

where y ~ \og{zi^/s/mj/^) is the rapidity of the vector meson, and <fi is the angle between 
p x and p 2 - Notice that the interference between the two mechanisms 7IP and IP7 is propor- 
tional to e\e 2 (p l ■ p 2 ) and introduces a charge asymmetry as well as an angular correlation 
between the outgoing protons. Clearly, the interference cancels out after integrating over </>, 
and the so integrated distributions will coincide for pp and pp collisions. 



B. Absorptive corrections 

We still need to correct for a major omission in our description of the production ampli- 
tude. Consider for example a rest frame of the proton 2 (the target) of the left panel of Fig[lJ 
Here, the virtual photon may be viewed as a parton of proton 1 (the beam), separated from 



1 Of course here the notation Q\ — ti applies only to the photon lines. 

2 While a trend towards s-channel helicity violating effects may be visible in the HI data [2j, they are surely 
negligible for our purpose, and within error bars, consistent with 1] and s-channel helicity conservation. 



4 



it by a large distance in impact parameter space. It splits into its cc-Fock component at a 
large longitudinal distance before the target, and to obtain the sought for production ampli- 
tude we project the elastically scattered cc system onto the desired J/\I/-final state (l5| . We 
entirely neglected the possibility [l6| that the photon's spectator partons might participate 
in the interaction, and destroy the rapidity gap(s) in the final state. Stated differently, for 
the diffractive final state of interest, spectator interactions do not cancel, and will affect the 
cross section. As a QCD-mechanism consider the interaction of a {cc}i{ggg}i-beam system 
with the target by multiple gluon exchanges (see Fig [2]). Then, for the J/\l/ final state 




FIG. 2: Left: the QCD two— gluon exchange mechanism for the Born-level amplitude. Right: a 
possible multigluon-exchange contribution that involves uncancelled spectator interactions. The 
impact parameters relevant for the discussion are indicated. 



of interest, the interaction of the cc-color singlet state is dominated by small dipole sizes 
r s ~ 4/mj/q, (the scanning radius of 13]). It can be exhausted by the minimal two-gluon 



color-singlet exchange, and will be quantifie d by the color-dipole cross section a(r) [15|, |18 
respectively its non-forward generalization (l9l 2p| . Let by be the tranverse separation of 
the J/\l/ and the target, and r the size of the cc-dipole as shown in Fig[2j Then, the 2^3 
amplitude of section III Al will involve, besides the vertex for the p — > 7p transition, the 
expectation value {J/^/\T^(r, by)]^) of 

r(°\r,b v ) = ±a(r)t N (b v ,B), (2.5) 

where t N (b,B) = exp(-6 2 /2B)/(2nB) is an optical density of the target. Systematic ac- 
count for the spectator interactions in QCD however is a difficult problem, as one cannot 
rely on the Abramovsky-Gribov-Kancheli (AGK) [2l[ cutting rules, when due account for 
color is taken j22j. To obtain at least a qualitative account of absorptive corrections we 
restrict ourselves to only a subclass of absorptive corrections, the 'diffractive cut', which 
contribution is model independent [23[. Regarding our {cc}i{ggg}i-system, with 6j denot- 
ing the constituent quarks' impact parameters, and b the impact parameter of the beam 
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proton, the absorbed amplitude in impact parameter space will contain 



r(r, b v , b) = - a(r) t N (b v , B) - -a(r)a qqq ({bi})t N (b v , B)t N {b, B el ) 

= rW(r, b v ) (l - l* qgq {{bi})t N (b, - r(°)(r, 6y) ■ 5 rf (6) . (2.6) 

In effect, we merely multiply the Born-level amplitude 1^°) by the probability amplitude for 
beam and target to pass through each other without inelastic interaction. In momentum 
space, we obtain the absorbed amplitude as 

M( Pl ,p 2 ) = J -0^S el (k)M(°\ Pl -k,p 2 + k) 

= M®(p 1 ,p 2 )-6M(p 1 ,p 2 ), (2.7) 

and with 

S el (k) = (2n)W 2 \k) - l -T{k) , T(k) = aZ(s) exp ( - l -B el k 2 ) , (2.8) 

the absorptive correction 5M reads 3 

5M( Pl ,p 2 ) = J ^- 2 T(k)M^( Pl -k,p 2 + k). (2.9) 

A number of improvements on this result can be expected to be relevant. Firstly, a more 
consistent microscopic treatment of spectator interactions along the lines of |22J would be 
desirable. Experience from hadronic phenomenology [25[ suggests that at Tevatron en- 
ergies, the purely elastic rescattering taken into account by eq. (12.91) are insufficient, and 
inelastic screening corrections will to a crude estimate lead to an enhancement of ab- 



sorptive corrections by a factor A ~ (a e i + (JD)/&ei [26J. Here an = Igsd + cdd, and 
&SD = cr(pp — > pX) , odd = cr(pp — > XY) are the cross sections for single-, and double- 
diffractive processes, respectively. Secondly, also the jp — > J/^p production amplitude 
will be affected by unitarity corrections. For example, with increasing rapidity gap Ay 
between J/\& and the target, one should account for additional s-channel gluons, and 
for sufficiently dense multiparton systems, the two-gluon exchange approximation for the 
7 — > J/ty transition used above, ultimately becomes inadequate. For relevant discussions of 
unitarity/saturation-effects in diffractive J/\l/-production, see j27[, the scaling properties of 
vector meson production in the presence of a large saturation scale are found in [20] . In our 



present approach, where the production amplitude is taken essentially from experiment, one 
must content oneself with the fact, that (some) saturation effects are effectively contained 
in our parametrization, and any extrapolation beyond the energy domain covered by data 
must be taken with great caution. 



III. RESULTS 

In this section we shall present results of differential cross sections for J/^> production. We 
shall concentrate on the Tevatron energy W = 1960GeV, where such a measurement might 



In the practical calculations below, for Tevatron energies, we take af^ t = 76 mb, B el = 17GeV -2 [2J 
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FIG. 3: A sketch of the elastic rescattering amplitudes effectively taken into account by eq.(J2 



be possible even at present. While in this paper we concentrate on the fully exclusive process 
pp — > ppJ/ty, pp — > ppJ/ty, it is important to realize, that from an experimental point of 
view there are additional contributions related to the exclusive production of Xc mesons 
and their subsequent radiative decays to J/^ffj. It may be difficult to measure/resolve the 
soft decay photons and therefore experimentally this contribution may be seen as exclusive 
production of J/if>. We note in this context that besides the scalar Xc(0 ++ ) meson, which 
exclusive production has been discussed in the literature (e.g. p and references therein), the 
axial-vector and tensor states Xc(l ++ ) and x c (2 ++ ) have larger branching fractions into the 
relevant channel. Although their exclusive production cross sections can be expected 



to be suppressed at low transverse momenta [g, [28|, a more detailed numerical analysis is 



not existing in the literature and would clearly go beyond the scope of the present paper. 



A. Distributions of J/i/j 

Let us start from the rapidity distribution of the J/^f shown in FigJH In the figure we 
present also the subsystem energies y/sijy/s^. At \y\ > 3 the energies of the 7p — ► J/ipp 
or 7p — > J/ipp subprocesses exceed the energy range explored at HERA. This may open 
a possibility to study J/ty photoproduction at Tevatron. This is interesting by itself and 
requires further detailed studies. In turn, this means that our estimate of the cross section 
far from midrapidity region requires extrapolations above the measured energy domain. In 
Fig[5]we collect rapidity distributions for different energies relevant for RHIC, Tevatron and 
LHC. We observe an occurence of a small dip in the distribution at midrapidities at LHC 
energy. The shape of the rapidity distribution at LHC energies however relies precisely on 
the above mentioned extrapolation of the parametrization of HERA-data to higher energies. 
Clearly a real experiment at Tevatron and LHC would help to constrain cross sections for 
•yp — > J I ipp process. 

In order to understand the origin of the small dip at midrapidity at LHC energy in 
FigE]we show separately the contributions of the two components (7lP,IP7 exchange) for 
Tevatron (left) and LHC (right). We see that at LHC energy the two components become 
better separated in rapidity. This reflects the strong rise of the J/^f photoproduction cross 
section with energy, which can be expected to slow down with increasing energy. Notice that 
the beam hadron hi moves along positive rapidities, so that, for example for the mechanism 
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FIG. 4: da/dy as a function of the J rapidity (y) for W = 1960 GeV. For a better understanding 
of the results we also show (dashed lines) the subsystem energies W\y = y/si and Wzv = \/~S2 in 
GeV. 

7IP it is the Pomeron exchange which 'propagates' over the larger distance in rapidity space. 
It would be interesting to confront our present simple predictions with the predictions of 
the approach which uses unintegrated gluon distributions - objects which are/were tested 
in other high-energy processes. This will be a subject of our forthcoming studies. 

Up to now we have not taken into account any restrictions on t\ and/or t 2 . In practice, 
it can be necessary to impose upper cuts on the transferred momenta squared. It is also 
interesting in the context of searches for odderon, to see how quickly the cross section for 
the "background" drops with t\ and t 2 . In Figj7]we show distribution in J/ip rapidity for 
different cuts on t\ and t 2 . Clearly, imposing a cut on t\ and t 2 removes the photon-pole 
contribution dominant at small momentum transfers. Even relatively small cut lowers the 
cross section considerably, and the dropping of the cross section is much faster than for 
the pomeron-odderon exchanges Imposing upper cuts on t\ and t 2 will therefore help 
considerably to obtain a possible "odderon-enriched" sample. We wish to repeat here that 
without absorption effects the rapidity distribution of J/ip in proton-proton and proton- 
antiproton collisions are identical. 

dcr(pp — > ppj/ty, W) da{pp — > ppj/^) 
dy dy 

It is interesting to stress in this context that it is not the case for transverse momentum 
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FIG. 6: da/dy as a function of y for the Tevatron and LHC energy. Individual processes are shown 
separately. Notice that the beam hadron hi of Fig. Q] moves at positive rapidities. No absorption 
corrections were included here. 



9 




FIG. 7: da/dy as a function of y for the Tevatron energy and different upper cuts on t\ and t^' 
t cut = 0.0 GeV 2 (solid), t cut = -0.05 GeV 2 (thin solid), t cut = -0.1 GeV 2 (dash-dotted), and 
t cut = -0.2 GeV 2 (dashed). 



distribution of J/ip, where 

dajpp -> ppj/m, W) dajpp -> ppj/^, W) 

This is demonstrated in FigJHl where we see, that at small transverse momenta of the vector 
meson, the interference enhances the cross section in pp collisions and depletes it in pp 
collisions. It is a distinctive feature of the mechanism discussed here, that vector mesons are 
produced with very small transverse momenta. The difference between proton-antiproton 
and proton-proton collisions survives even at large rapidities of J/ip. When integrated over 
the J / ip transverse momentum, and in absence of absorptive corrections, cross sections will 
again be identical in the pp and pp cases. 



B. Distributions of (anti— ) protons 

Now we shall proceed to distributions related to (anti-) protons. In Figj9]we show distri- 
butions in the transferred momenta squared (identical for t\ and t 2 ). We show separately the 
contributions of 7IP and IP7 exchanges. The figures clearly display the strong photon-pole 
enhancement at very small t. 

In order to better understand the distributions in t\ or ti in Fig JlO] we show how t\ and ti 
are correlated. Here we do not make any restrictions on the rapidity range. The significant 
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FIG. 8: The distribution da/dp^ of J/ip MS cX function of J ftp transverse momentum for different 
intervals of rapidity: —0.5 < y < 0.5 (left panel) and 3.5 < y < 4.5 (right panel) at W = 1960 
GeV. The result for pp collisions is shown by the solid line and the result for pp collisions by the 
dashed line. No absorption corrections were included here. 




FIG. 9: da/dti/2 as a function of Feynman t l / 2 for W = 1960 GeV. The photon-exchange (dashed) 
and pomeron-exhange (dotted) contributions are shown in addition. No absorption corrections were 
included. 
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enhancements of the cross section in the form of ridges along t± ~ and t 2 ~ are again due 
to the massless photon-exchange, and most of the integrated cross section comes from these 
regions. The pomeron-odderon and odderon-pomeron exchange contributions considered in 
Ref. [12[| would not exhibit such a significant local enhancements and would be smeared over 
broader range in the (£1,^2) space. Therefore in the dedicated searches for the odderon 
exchange upper cuts on t 1 and t 2 should be imposed, and t upper = —0.2 GeV seems to be a 
good choice. 




FIG. 10: Two-dimensional distribution in t\ and £2 for the Tevatron energy W = 1960 GeV. In 
this calculation a full range of the J/tp rapidities was included. No absorption corrections were 
included here. 

We repeat that the reaction considered in this paper leads to azimuthal correlations 
between outgoing proton and antiproton. In FigJTTJ we show the corresponding angular 
distribution for proton-antiproton collision (solid line). For reference we also show, by the 
dotted line, the incoherent sum of the 7IP and IP7 mechanisms. The distribution for proton- 
proton collisions (dashed line) is shown for comparison alsoat the Tevatron energy. Clearly 
the interference terms in both reactions are in opposite phase due to different electric charges 
of proton and antiproton. In the absence of absorptive corrections, we have 

^ = A±Bcos<j> (3.3) 
d(j) 

for pp (+) and pp (-) collisions respectively. The interference effect (B/A) here is at the level 
of ~ 40-50 %. 

In Fig{12] we show the two-dimensional distributions differentially in both rapidity and 
azimuthal angle. Interestingly, the interference effect is significant over broad range of 
J/ip rapidity, which is reflected in the fact that even at large J/ip rapidities one observes 
ansisotropic distributions in the azimuthal angle. 

Up to now we have considered only spin-preserving contributions. Now we wish to show 
the effect of electromagnetic spin-flip discussed in the previous section. In Fig{13] we show 
the ratio of helicity-flip to helicity-preserving contribution. The ratio is a rather flat function 
of ti and t 2 . At t\ = —1 GeV 2 and t 2 = —1 GeV 2 the ratio reaches about 0.4. 
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FIG. 11: da/dtp as a function of (f> for W = 1960 GeV. The solid line corresponds to a coherent 
sum of amplitudes whereas the dashed line to incoherent sum of both processes. No absorption 
corrections were included here. 




FIG. 12: do /dydcj) for W = 1960 GeV and for pp (left panel) and pp (right panel) collisions. No 
absorption corrections were included here. 
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FIG. 13: The ratio of helicity-flip to helicity-preserving contribution as a function of t\ and t%. 
C. Absorption effects 

Now we will show the effect of absorptive corrections discussed in section III Bl on various 
differential distributions. 

Let us start from the presentation of the effects of absorption for selected points in phase 
space. In FigJTHwe show the fully differential cross section da/dydtidt 2 d<j) as a function of <f) 
for selected (fixed) values of ti, t 2 and for y = 0. We show results for pp (left panel) and pp 
(right panel) collisions for the same center-of-mass energy W = 1960 GeV. While at smaller 
ti t 2 we observe a smooth reduction of the Born-level result, absorptive corrections induce 
a strong ^-dependence at larger ti i2 - The positions of the diffractive minima which appear 
as a consequence of cancellations of the Born and rescattering amplitudes move with the 
value of t = ti — t 2 . In Fig{T5]we present fully differential cross section as a function of t\ 
and t 2 for y=0 and <fi = ir/2. For proton-proton scattering we observe clearly a diffractive 
minimum for ti — t 2 ~ 0.4 GeV 2 . In Fig. [in] we show the fully differential cross section as 
a function of the transverse momentum squared p 2 at a fixed value p 2 = lGeV 2 at rapidity 
y = 0. The rich structure as a function of transverse momenta and azimuthal angles is 
also revealed by this plot. The plots in FigJT7]give an idea, to which extent the diffractive 
dip-bump structure depends on the details of our treatment of absorption. Here we show, 
by the dotted line, the cross section calculated for the Born-level amplitude. The solid line 
shows the result with elastic scattering included, and for the dashed and dash-dotted lines 
we enhanced the rescattering amplitude T of section HlBl by a factor A = 1.2 and A = 1.5 
respectively. The region of very small p 2 is entirely insensitive to rescattering, reflecting the 
ultraperipheral nature of photon exchange. The diffractive dip-bump structure, situated at 
larger transverse momenta reveals a dependence on the strength of absorptive corrections. 
This concerns the position of dips as well as the strength of the cross section in various 
windows of phase space. 

The sensitivity to rescattering is however washed out in integrated observables - clearly 
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FIG. 14: Fully differential cross section da / ' dydt\dtid<\) as a function of <p for y=0 and different 
combinations of t\ = t2 (—0.1, —0.3, —0.5 GeV 2 (from top to bottom)) for pp (left panel) and pp 
(right panel) reactions. The solid lines include rescattering, while the dashed lines correspond the 
Born-level mechanism only. 




FIG. 15: The fully differential cross section da / 1 dydt\dtid<\) as a function of t\ and t-x for y=0 and 
4> = vr/2 for pp (left panel) and pp (right panel) reactions. 

the contribution from low transverse momenta, which is not strongly affected by rescattering, 
is large. This becomes apparent in Figs. [13, [19] and [20j 

In Fig[l8]we show the ratio of the cross section with absorption to that without absorption 
as a function of t\ and t2, again for pp (left) and pp (right). As a consequence of averaging 
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FIG. 16: Fully differential cross section da j 1 dydp\dp\ d(p as a function of p\ at y = and 
p\ = lGeV 2 for pp (left panel) and pp (right panel) collisions at W = 1960 GeV. Absorptive 
corrections (elastic rescattering) are included. Solid, dashed, and dash-dotted lines refer to different 
values of the azimuthal angle (j) between outgoing (anti-)protons. 



over different phase-space configurations the diffractive minima disappeared. Again the 
plot reflects, that due to the ultraperipheral nature of the photon-exchange mechanism, 
absorption is negligible at very small tip, and rises with ti and/or t 2 - On average, absorption 
for the pp reaction is smaller than for the pp case. 

It is important to stress again, that the absorptive corrections in differential cross sections 
cannot be accounted for by simply a constant suppression factor, but show a lively depen- 
dence over phase space. In Fig. [TH]we show the differential cross section da/dt. The dashed 
line shows the result without absorption, the solid lines include absorptive corrections. They 
differ by a factor 1.5, by which rescattering had been enhanced in the lower curve. This 
enhancement of rescattering shows a modest effect, quite in agreement with the expectation 
mentioned above. The dependence of absorption on t is quantified by the ratio of full and 
Born-level cross sections shown by the dotted lines. 
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FIG. 17: Fully differential cross section da / dydpfdp^dcp as a function of p\ at y = and p\ = 
lGeV 2 for pp (left panel) and pp (right panel) collisions at W = 1960 GeV. The azimuthal angle 
(j) is taken = in the top row and (f> = tt in the second row. Shown by the dotted curve is the 
Born-level cross section, without absorptive corrections included. The solid curve shows the result 
with elastic rescattering between (anti-) protons included, while rescattering has been enhanced by 
a factor A = 1.2 in the amplitude for the dashed curve, and by a factor A = 1.5 for the dash-dotted 
curve. 



In Fig. [20] we show the suppression factor 

< b {Pv} >_ da^/dpldy ' {6A} 

as a function of the J/\l/ transverse momentum at y = 0. It is important to emphasize once 
more the strong functional dependence on p v , which is different for pp and pp collisions. 
Again, playing with the strength of absorptions shows a modest effect. 

The different behaviour of absorptive corrections in pp and pp collisions is an interesting 
observation. It derives from the fact that rescattering corrections lift the cancellation of 
the interference term after azimuthal integration. Finally, let us comment on the expected 
reduction of rapidity distributions from absorptive corrections. These are, finally, rather flat 
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FIG. 18: The ratio of the cross sections with absorption to that without absorption for pp (left 
panel) and pp (right panel) scattering. Here the integration over y and 4> was performed. 




FIG. 19: Differential cross section da/dt, integrated over — 1 < y < 1 for pp collisions. Dashed line: 
no absorptive corrections. Solid lines show the cross section with absorptive corrections included. 
Upper solid line: purely elastic rescattering. Lower solid line: rescattering enhanced by a factor 
A = 1.5 in the amplitude. The dotted lines show the respective ratios d a Born + Rescatt - /da Born . 



functions of y. For the ratio 



we obtain, in pp collisions 



.,, rf(T Born+Rescatty dy 



< S\y) >= 



da Born /dy 



<S 2 (y = 0)> «0.9, < S 2 (y = 3) > 

pp 



0.8, 



pp 



(3.5) 



(3.6) 
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and for pp collisions 



< S 2 {y = 0) > « 0.85 , < S 2 (y = 3) > 



0.75. 



We note that this results in a small charge asymmetry 

da(pp) /dy — da{pp) / dy 
da (pp) /dy + da (pp) / dy 

which derives entirely from absorptive corrections. 



2-3% 



(3.7) 



(3.8) 



IV. CONCLUSIONS 

In this paper we have calculated differential cross sections for exclusive J/tp production 
via photon-pomeron(7lP) and pomeron-photon (IP7) exchanges at RHIC, Tevatron and 
LHC energies. Measurable cross sections were obtained in all cases. We have obtained an 
interesting azimuthal-angle correlation pattern due to the interference of the " TP and TP 7 
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mechanisms. The interference effect survives over almost the whole range of J/ip rapidities. 
At the Tevatron energy one can potentially study the exclusive production of J/ip at the 
photon-proton center-of-mass energies 70 GeV < W 1P < 1500 GeV, i.e. in the unmeasured 
region of energies, much larger than at HERA. At LHC this would be correspondingly 200 
GeV < W 1P < 8000 GeV. At very forward rapidities this is an order of magnitude more than 
possible with presently available machines. Due to the photon-pole, the differential cross 
section is concentrated in the region of very small ti or/and ti- Imposing cuts on t\ and 
t2 lowers the cross section considerably. Electromagnetic helicity-flip processes play some 
role only when both \t\\ and ^1 are large, that is in a region where also the hypothetical 
hadronic, Odderon exchange, contribution can be present. It is a distinctive feature of the 
production mechanism, that mesons are produced at very small transverse momenta, where 
the interference of 7IP and IP7 mechanisms induces a strongly different shape of vector- 
meson py- distributions in pp vs. pp collisions. We also estimated absorption effects on 
various distributions. In some selected configurations the absorptive corrections lead to the 
occurence of diffractive minima. Naturally, the exact place of diffractive minima depends 
on the values of the model parameters, but they are washed out when integrated over the 
phase space or even its part. Absorptive corrections for differential distributions are lively 
functions of transverse momenta, and cannot be accounted for simply by constant suppres- 
sion factors. We have found that absorptive corrections induce a small charge-asymmetry 
in rapidity distributions and total production cross sections. In the present paper we have 
concentrated on exclusive production of J / ip at energies ^ > 200 GeV. The formalism used 
here can be equally well applied to exclusive production of other vector mesons, such as (ft, 
T as well as to the lower energies of e.g. FAIR, J- PARC, RHIC. For J/^-production, espe- 
cially recent parametriztions [29( of the photoproduction cross section from threshold to the 
highest energies may prove useful. We leave such detailed analyses for separate studies. The 
processes considered here are also interesting in the context of recently proposed searches 
for identifying the Odderon. We find that the region of midrapidities (— 1 < y < 1) and 
ti,t2 < —0.2 GeV 2 seems the best in searches for the odderon exchange. Should data reveal 
deviations from the conservative predictions given by us, a detailed differential analysis of 
both, photon and odderon-echange processes, including their interference, in y,t\,t2,4> will 
be called upon. 
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VI. APPENDIX 

In Ref. [2[ the differential cross section ^ for the reaction 7*p — > J / ipp was parametrized 

as 

>,,,«■)=!! (wy™-\ PiBot) (^* v (6 . 1} 

at dt t=o,w=w \W J \ m j/ip + Q I 
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where a(t) = ao + a!t . The values of parameters found from the fit to the data are: 
^\t=o,w=w = 326 nb/GeV 2 , W = 95 GeV, B = 4.63 GeV" 2 , a = 1.224, a' = 0.164 
GeV"' 2 , n = 2.486. 

Assuming the dominance of the helicity-conserving transitions, and neglecting the real 
part, one can write 

I j / \ a(t)-i / m 2 \ n ^ 2 

M(s,t,Q 2 ) = 5 XjXv 5 XpX pl is\ /167T-J exp( J B t/2) 2 J ^ n2 

(6.2) 

identical for each combination of particle helicities. In our case of hadropro duct ion the 
amplitude is a function of either (si,ii,Q|) or (s 2 ,i2,Qi)- 
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